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I.  INTRODUCTION 

Recent  involvement  in  the  SENTRY  program  and  expected  future  involvement 
in  the  Advanced  Tactical  Missile  (ATM)  program  prompts  this  technical 
assessment  of  the  application  of  Richman's[l]  Q-function  to  waveform  design. 
Because  of  its  nature,  the  Q-function  has  generally  been  limited  to  waveform 
design  in  the  Ballistic  Missile  Defense  (BMD)  community.  However,  even  in 
this  limited  application,  it  has  not  always  been  correctly  applied.  Because 
the  Q-function  derives  from  Woodward ' s [ 2 ]  ambiguity  function,  the  first  part  of 
this  work  will  address  the  ambiguity  function  and  point  out  some  misconcep¬ 
tions  concerning  this  function.  Following  this  brief  diversion  a  general 
assessment  of  the  Q-functlon  will  be  made. 

II.  WOODWARD'S  AMBIGUITY  FUNCTION 

Generally,  the  purpose  of  a  radar  transmission  is  to  detect  an  object  and, 
in  many  cases,  to  furnish  other  information  concerning  that  object.  For  a 
given  radar  dwell,  this  object,  or  radar  target,  is  .located  in  range-doppler 
or,  equivalently,  T-f^  space.  (Other  location  information  such  as  angle  may 
be  available,  but  this  Information  is  not  waveform  dependent).  The  basic 
question  considered  by  Woodward  was  the  resolution  of  two  targets  residing  at 
different  points  in  r-fj  space.  Then  the  problem  was  to  design  the  wavfeform 
so  as  to  maximize  the  difference  in  some  sense  between  the  returns  of  the  two 
distinct  targets. 

r 

Therefore,  consider  a  point-target  located  at  the  origin  (0,  0)  of 
r-fd  space  and  a  second  point-target  located  arbitrarily  at  (T,f<j)  in  this 
space.  The  waveform  received  from  the  target  located  at  the  origin  is  iden¬ 
tical  to 'the  transmitted  complex  waveform 

j2itf0t  . 

Z(t)  -  u(t)eJ  °  (1) 


The  time-delayed  and  Doppler-shifted  waveform  received  from  the  second 
target  is 

j2*(f0-fl»)(t-c) 

Zj(t)  -  u(t--t)e  .  (2) 

A  convenient  measure  of  the  difference,  or  distance,  between  these  two  wave¬ 
forms  is  the  integrated  square  magnitude 


E2  -/|  *(t)-*1(t)  |  2dt.  (3) 

Then  the  waveform  which  best  resolves  the  two  targets  is  the  one  which 
maximizes  E2.  Rewriting: 

e2  -  /|„<t).32,fot  -  |  2dt 

-  /[u(t)u*(t>  -  <i«<t)u<t-T)e"J2*<f°'fd)VJ2”fdt 

-  u(t)»*<t-,>,J2’<'°'fd)V2,,f',t  ♦ 

■  /[u(t)u*(t)  +  u(t-T)u*(t-t)ldt 


-/I  u(t)u*(t-x) 


j2w<f0-fd)x  j 2 xtf 
e  e  dt 


+ 


(u(t)u*(t-x)e 


j2n(f0-fd)t  j2nfdt 
e 


* 

)  1  dt 


-2/|  u(t)  |  2dt  -  2Re  jeJ2ll(f°“fd)*  Ju(t)u*(t-x)  eJ2*fdtdt)  (4) 

where  *  signifies  complex  conjugate. 

The  first  term  in  Eq.  (4)  is  just  twice  the  energy  of  the  transmitted 
waveform  and  will  normally  be  maximized  for  the  particular  application  under 
investigation.  Therefore,  in  order  to  maximize  E?  it  is  necessary  to  minimize 
the  magnitude  of  the  last  term  of  Eq.  (4).  Thus,  it  is  desired  to  minimize 

|  |  "  |Ju(t)u*(t-x)e^2  d  dt  |  (5) 

or,  equivalently,  to  minimize 

|  x(*»fd)  |  2  *  |  Ju(t)u*(t-x)e^  d  dt  |  2  (6) 

which  is  Woodward's  ambiguity  function.  It  is  common  practice  to  normalize 
Eq.  (5)  or  (6)  so  that 

|  X<i,fd)  |  /  |  X<0,  0)  |  (7) 


or 


|  x(x,fd)  I  2/  I  x<0,0)  |  2  .  T(x,fd)  (8) 

are  the  functions  with  which  one  normally  works. 

There  are  two  Important  points  to  note  here.  First,  one  observes  that 
x(x,fd)  coincidentally  has  the  form  of  an  ideal  matched  filter  response. 

Thus,  the  ambiguity  function  has  often,  but  erroneously,  been  referred  to  as 
the  matched  filter  response  of  a  system.  Starting  from  this  point,  research¬ 
ers  have  constructed  mismatched  receiver  responses  and  called  these  responses 
cross-ambiguity  functions.  The  second  point  is  that  the  ambiguity  function  is 
totally  processor  independent.  Therefore,  a  cross-ambiguity  function  has  no 
meaning  whatsoever  in  the  Woodward  sense.  That  is,  the  radar  processor  will 
in  no  way  affect  the  waveform  ambiguity  function. 

III.  RICHMAN'S  Q- FUNCTION 

Rlchman  Introduced  a  function  defined  by 


Q<fd> 


T(x,fd)dx 


(9) 


which  has  commonly  come  to  be  called  a  Q-function 


More  practically,  one  runs  into  a  relative  Q-function  defined  by 


Q(fd> 


fd)dT 

IV^T 


(10) 


where  is  the  ambiguity  function  of  a  single  pulse  of  a  burst  waveform 

with  ambiguity  function  ?(T,fd)<  Since  the  same  notation  Q(f<j)  is  used  for 
both  of  these  forms,  one  should  be  careful  in  using  functional  data  supplied 
by  someone  else* 


An  integrated,  properly  weighted  ambiguity  function  results  in  a  quantity 
proportional  to  clutter  power.  Therefore,  the  relative  Q-function  is  a 
measure  of  the  improvement  in  clutter  rejection  of  a  burst  waveform  over  that 
provided  by  a  single  pulse^  of  the  burst  when  the  clutter  is  uniformly  extend¬ 
ed  in  range.  This  is  precisely  the  type  of  clutter  exhibited  by  a  wide-spread 
rain  storm.  In  a^WND  environment,  the  volume  filled  by  tank,  breakup  clutter 

closely  approximates  extended  uniform  clutter. 

♦ 

Justification  for  utilization  of  the  Q-function  is  straightforward.  If 
operation  is  in  range-extended  clutter,  then  target  discrimination  must  be 
based  upon  Doppler  information.  Therefore,  the  waveform  ambiguity  function 
should  exhibit  small  amplitudes  at  Doppler  values  away , from  the  central  peak. 
Additionally,  if  the  ran^e  extended  clutter  is  essentially  homogeneous,  then  a 
one  dimensional  measure  of  the  waveform's  effectiveness  may  be  obtained  by 
collapsing  the  ambiguity  function  in  range.  It  can  be  seen  that  the  function 
given  by  Eq.  (10)  account's  both  for  the  increased  clutter  power  resulting  from 
the  transmission  of  a  burst  waveform  as  well  as  the  decrease  in  clutter  power 
resulting  from  Doppler  filtering. 


IV.  UTILIZATION  OF  THE  Q-FUNCTION 


In  an  application  such  as  BMD,  where  the  search  and  track  functions  are 
separated,  a  great  deal  of  flexibility  and  optimization  may  be  achieved  by 
having  a  variety  of  waveforms  available  for  transmission.  On  the  other  hand, 
a  track-while-scan  system  has  limited  flexibility  in  regard  to  waveform  design 
since  both  detection  and  parameter  estimation  generally  must  be  accomplished 
simultaneously. 

Equation  (10)  gives  a  measure  of  clutter  rejection  provided  by  a  burst 
waveform  over  a  single  pulse.  This  clutter  rejection  is  based  solely  upon 
Doppler  discrimination.  A  uniform  amplitude  burst  will  give  some  clutter 
rejection  Improvement  over  a  single  pulse,  but  a  weighted  burst,  because  of 
its  reduced  Doppler  sidelobes,  will  produce  a  much  more  dramatic  clutter 
rejection.  Unfortunately,  this  weighting  comes  about  at  the  expense  of 
transmitted  power;  however,  some  clutter  environments  are  so  severe  that  wave¬ 
form  amplitude  tailoring  is  a  necessity. 

The  Q-function  does  have  some  utility  in  evaluating  various  weighted 
waveforms;  however,  this  application  is  somewhat  limited  since  an  exotic 
weighting  scheme  will  produce  little  improvement  over  standard  weighting 

Ilf  the  individual  pulses  have  different  bandwidths,  then  the  pulse  with  maxi¬ 
mum  bandwidth  is  chosen  for  comparison. 


schemes  such  as  Taylor's.  Some  clutter  rejection  Improvement  Is  available  by 
varying  the  pulse  widths  In  a  burst.  This  Improvement  appears  to  be  much  less 
than  that  obtainable  by  weighting  but  transmission  power  Is  conserved.  The 
Q-function  may  be  much  more  useful  in  evaluating  pulse  width  modulated  wave¬ 
forms  than  in  evaluating  weighted  waveforms.  However,  pulse  width  modulation 
produces  many  processing  problems  and  is  seldom,  if  ever,  employed. 

V.  CONCLUSION 

The  Q-function  has  limited  application  and  is  useful  primarily  in  eval¬ 
uating  time-domain  weighting  performance.  It  is  seldom  used  to  compare  two 
different  waveforms;  that  is,  where  the  individual  pulses  are  different  for 
the  two  waveforms.  For  the  most  part,  the  many  different  weighting  schemes 
have  been  over-evaluated  many  years  ago  and  should  now  be  selectable  via  engi¬ 
neering  judgment  rather  than  complicated  analysis.  The  structure  of  the 
single  pulse  from  which  the  burst  is  constructed  appears  to  have  little  effect 
on  the  Q-function  and  should  be  selected  on  other  basic  considerations.  Out¬ 
side  the  BMD  community,  at  least  one  weighted  waveform  should  be  strongly  con¬ 
sidered  for  use  during  periods  of  heavy  weather  activity.  The  use  of  more 
than  one  special  waveform  should  be  strongly  justified. 
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APPENDIX  A 


Computer  Generation  of  Burst  Ambiguity  and  Q-Function 

The  ambiguity  and  Q-functions  are  directly  obtainable  from  y(T,fd)  given 
by  Eq.  (5).  Let  u(t)  be  a  single,  unity  amplitude  pulse  of  a  burst  waveform 
consisting  of  N  such  pulses  with  scaled  amplitudes  and  pulse  repetition  period 
AT  as  shown: 


+•  AT 


In  the  following  a  receiver  response  function  will  be  calculated  mainly 
because  of  its  widespread  use  (also  misnomered  as  a  cross-ambiguity  function). 
The  actual  ambiguity  function  is  obtained  where  the  transmit  and  receive 
weights  are  a  complex  conjugation  of  each  other.  Therefore,  let 

N 

uL(t)  -  Z  Aku[t-(k-I)ATj 
k-1 
N 

u2(t)  -  I  B*u[t-(A-1)AT] 

A«1 

N 

u2*(t)  =  Z  B**u*[t-(A-1)AT]  (A-l) 

A- 1 

where  Ak  and  Bj*  are  the  transmit  and  receive  weights,  respectively. 

Then* 

r  j2nfdt 

y(t,f(j)  ui(t)u2*(t-x)e  dt 

fN  N  J2itfdt 

=  1  Z  Aku[t-(k-l)AT]  Z  Bj^*u*  [  t-(  A-l ) AT-t] e  dt 

A-l  A-l 

N  N  40  f  ► 

/1 2itfdt 

u( t-(k-l ) AT) u*l t-( A-l )AT-t] e  dt  (A-2) 

A.- I 

^Here  the  nomenclature  for  the  ambiguity  function  is  used  although  the  result 
is  an  ambiguity  function  only  where  transmit  and  receive  weighting  are  complex 
conjugates. 
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Letting  x  -  t-(k-l)AT 


N  N 

X(x.fd)  “  s  £  AkBj*  Ju(x)u*(x-{(  *.-k)AT+x}  ]  e 


j2itfdx  j2itfd(k-l)AT 

e  dx 


k-1  A-l 
N 


N 


j2itfd<k-l)AT 

-  E  Ak  e  E  Bi*X1[(i-k)AT  +  t,fd] 

k=l  X-l 

where  X^(<)  Is  the  single  pulse  ambiguity  function. 

Discretizing  x  and  fd  results  ini 


(A-3) 


x(mAT.nAfd) 


k-1  *-l 


+  mAx.nAfdl,  m,n  =  0,1,... 


(  A-4) 


Thus  the  receiver  response  is  just  a  linear  combination  of  shifted  single 
pulse  ambiguity  functions.  It  remains  only  to  calculate  this  single  pulse 
ambiguity  function: 

Xi(x,fd)  -  u(t)u*(t— c)e^  ^  dt  (A”5) 
Discretizing  and  approximating  the  integral  by  a  summation: 


N-l 

Xj(mAt,nAfd)  *»  I 

i-0 


e 


jn!2it 

N 


(A-6) 


Omitting  the  scale  factor,  the  result  is  just  an  inverse  discrete  Fourier 
transform  for  each  m  and  is  easily  calculated  using  fast  Fourier  transform 
(FFT)  techniques.  Therefore 


X^mAr.kAfd)  -■  IFFT 


“*(i5E  '  “T)] 


(A-7) 


The  procedure  for  calculating  Eq.  (A-3)  is: 

1.  Select  an  appropriate  sampling  interval  for  the  single  pulse  (Ex.  O.lps). 

2.  Select  the  desired  number  of  Doppler  samples  (Ex.  4096)2. 

3.  Select  a  value  of  x  equal  to  some  multiple  of  the  sampling  interval. 

4.  Form  the  product  array  u(i)  u*  (i-x)  consisting  of  Nj  nonzero  samples. 


^Note  that  yjl(  i-k)  AT+mAx,nAf(j]  -  0  for  -AT>l(S.-k)AT+mAx]>AT 

o  T  T 

‘This  determines  Tp  in  Figure  A-2  since  _R  -  sampling  interval  or  .  P  ■  0.1 

giving  Tp  •  409.6  ps  here.  N  4096 


5.  Append  zeros  to  this  array  to  total  N  (4096)  samples. 

6.  Compute  the  inverse  FFT  of  this  array. 

7.  Result  is  doppler  values  (with  ~  Hz  spacing)  for  this  choice  of  x. 

TP 

8.  Repeat  for  all  desired  values  of  x. 

With  the  arrays  resulting  from  the  above  procedure,  the  receiver  response 
is  now  easily  calculated  using  Eq.  (A-4).  The  Q-function  is  obtainable  by 
simply  collapsing  the  ambiguity  function  in  range.  However,  do  not  forget  to 
apply  the  proper  normalization  to  both  the  ambiguity  and  Q-function. 

Figures  A-l  and  A-2  give  examples  of  a  Q-function^.  Figure  A-l  is  the 
Q-function  of  a  filled  uniformly  weighted  16-pulse  burst  where  each  pulse  is 
linearly  frequency  modulated  with  a  bandwidth  of  30  MHz.  Figure  A-2  is  iden¬ 
tical  to  A-l  except  that  a  50  dB  Chebyshev  weighting  has  been  applied  across 
the  burst.  These  figures  correspond  to  Eq.  (9)  except  that  a  different  nor¬ 
malization  than  that  given  by  Eq.  (8)  was  used. 


^Figures  supplied  by  Dr.  Jerry  Cross,  Teledyne  Brown  Engineering,  Huntsville, 
Alabama. 


•AN0MOTH  -  30.00 MHZ  PULSC  LENGTH  •  I  SO MtCMOKCONM 


Figure  A-l.  Filled,  uniformly  weighted  16-pulse  burst  Q-function. 
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